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This application is also a continuation-in-part of and claims priority to U.S. patent 
application serial number 10/633,186 filed July 7, 2003 which claims priority to U.S. 
provisional patent application serial number 60/423,881 filed November 5, 2002 and U.S. 
provisional patent application serial number 60/423,881 filed August 2, 2002, and to U.S. 
5 patent application serial number 10/443,126 filed May 20, 2003 which claims priority to U.S. 
provisional patent application serial number 60/423,881 filed November 5, 2002 and to U.S. 
provisional application 60/400,624 filed August 2, 2002. 

This application also claims priority to U.S. provisional patent application serial 
number 60/470,525 filed May 12, 2003, and to U.S. patent application serial number 
10 10/165,556 filed June 7, 2002. All of the above applications are herein incorporated by 
reference in their entirety as if fully set forth herein. 

Field of the Invention 

( This invention pertains to the field of obstetrics, particularly to ultrasound-based non- 

1 5 invasive obstetric measurements. 

Background of the Invention 

Measurement of the amount of Amniotic Fluid (AF) volume is critical for assessing 
the kidney and lung function of a fetus and also for assessing the placental function of the 
20 mother. Amniotic fluid volume is also a key measure to diagnose conditions such as 
polyhydramnios (too much AF) and oligohydramnios (too little AF). Polyhydramnios and 
oligohydramnios are diagnosed in about 7-8% of all pregnancies and these conditions are of 
concern because they may lead to birth defects or to delivery complications. The amniotic 



25315 

PATENT TRADEMARK OFFICE 



-2- 



DXUC-I-I020AP 



Black Lowe & Graham f 

ait- 



816 Second Avenue 
Seattle, Washington 98104 
206.381.3300 • F: 206.381.3301 



fluid volume is also one of the important components of the fetal biophysical profile, a major 
indicator of fetal well-being. 

The currently practiced and accepted method of quantitatively estimating the AF 
volume is from two-dimensional (2D) ultrasound images. The most commonly used measure 
5 is known as the use of the amniotic fluid index (AFI). AFI is the sum of vertical lengths of 
the largest AF pockets in each of the 4 quadrants. The four quadrants are defined by the 
umbilicus (the navel) and the linea nigra (the vertical mid-line of the abdomen). The 
transducer head is placed on the maternal abdomen along the longitudinal axis with the 
patient in the supine position. This measure was first proposed by Phelan et al (Phelan JP, 
10 Smith CV, Broussard P, Small M., "Amniotic fluid volume assessment with the four- 
quadrant technique at 36-42 weeks 1 gestation," J Reprod Med Jul; 32(7): 540-2, 1987) and 
then recorded for a large normal population over time by Moore and Cayle (Moore TR, 
Cayle JE. "The amniotic fluid index in normal human pregnancy," Am J Obstet Gynecol 
May; 162(5): 1168-73, 1990). 

15 Even though the AFI measure is routinely used, studies have shown a very poor 

correlation of the AFI with the true AF volume (Sepulveda W, Flack NJ, Fisk NM., "Direct 
volume measurement at midtrimester amnioinfiision in relation to ultrasonographic indexes 
of amniotic fluid volume," Am J Obstet Gynecol Apr; 170(4): 1 160-3, 1994). The correlation 
coefficient was found to be as low as 0.55, even for experienced sonographers. The use of 

20 vertical diameter only and the use of only one pocket in each quadrant are two reasons why 
the AFI is not a very good measure of AF Volume (AFV). 

Some of the other methods that have been used to estimate AF volume include: 
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Dye dilution technique. This is an invasive method where a dye is injected into the 
AF during amniocentesis and the final concentration of dye is measured from a sample of AF 
removed after several minutes. This technique is the accepted gold standard for AF volume 
measurement; however, it is an invasive and cumbersome method and is not routinely used. 

5 Subjective interpretation from ultrasound images. This technique is obviously 

dependent on observer experience and has not been found to be very good or consistent at 
diagnosing oligo- or poly-hydramnios. 

Vertical length of the largest single cord-free pocket. This is an earlier variation of 
the AFI where the diameter of only one pocket is measured to estimate the AF volume. 

10 Two-diameter areas of the largest AF pockets in the four quadrants. This is 

similar to the AFI; however, in this case, two diameters are measured instead of only one for 
the largest pocket. This two diameter area has been recently shown to be better than AFI or 
the single pocket measurement in identifying oligohydramnios (Magann EF, Perry KG Jr, 
Chauhan SP, Anfanger PJ, Whitworth NS, Morrison JC, "The accuracy of ultrasound 

15 evaluation of amniotic fluid volume in singleton pregnancies: the effect of operator 
experience and ultrasound interpretative technique," J Clin Ultrasound, Jun; 25(5):249-53, 
1997). 

The measurement of various anatomical structures using computational constructs are 
described, for example, in U.S patent 6,346,124 to Geiser, et al. (Autonomous Boundary 
20 Detection System For Echocardiographic Images). Similarly, the measurement of bladder 
structures are covered in U.S. patent 6,213,949 to Ganguly, et al. (System For Estimating 
Bladder Volume) and U.S. patent to 5,235,985 to McMorrow, et al., (Automatic Bladder 
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Scanning Apparatus). The measurement of fetal head structures is described in U.S. patent 
5,605,155 to Chalana, et al., (Ultrasound System For Automatically Measuring Fetal Head 
Size). The measurement of fetal weight is described in U.S. patent 6,375,616 to Soferman, et 
al. (Automatic Fetal Weight Determination). 

5 Pertaining to ultrasound-based determination of amniotic fluid volumes, Segiv et al. (in 

Segiv C, Akselrod S, Tepper R., "Application of a semiautomatic boundary detection 
algorithm for the assessment of amniotic fluid quantity from ultrasound images." Ultrasound 
Med Biol, May; 25(4): 515-26, 1999) describe a method for amniotic fluid segmentation from 
2D images. However, the Segiv et al. method is interactive in nature and the identification of 
10 amniotic fluid volume is very observer dependent. Moreover, the system described is not a 
dedicated device for amniotic fluid volume assessment. 

Grover et al. (Grover J, Mentakis EA, Ross MG, "Three-dimensional method for 
determination of amniotic fluid volume in intrauterine pockets." Obstet Gynecol, Dec; 90(6): 
1007-10, 1997) describe the use of a urinary bladder volume instrument for amniotic fluid 

15 volume measurement. The Grover et al. method makes use of the bladder volume instrument 
without any modifications and uses shape and other anatomical assumptions specific to the 
bladder that do not generalize to amniotic fluid pockets. Amniotic fluid pockets having shapes 
not consistent with the Grover et al. bladder model introduces analytical errors. Moreover, the 
bladder volume instrument does not allow for the possibility of more than one amniotic fluid 

20 pocket in one image scan. Therefore, the amniotic fluid volume measurements made by the 
Grover et al. system may not be correct or accurate. 
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None of the currently used methods for AF volume estimation are ideal. Therefore, 
there is a need for better, non-invasive, and easier ways to accurately measure amniotic fluid 
volume. 

Summary of the Invention 

5 The preferred form of the invention is a three dimensional (3D) ultrasound-based 

system and method using a hand-held 3D ultrasound device to acquire at least one 3D data 
set of a uterus and having a plurality of automated processes optimized to robustly locate and 
measure the volume of amniotic fluid in the uterus without resorting to pre-conceived models 
of the shapes of amniotic fluid pockets in ultrasound images. The automated process uses a 
10 plurality of algorithms in a sequence that includes steps for image enhancement, 
segmentation, and polishing. 

A hand-held 3D ultrasound device is used to image the uterus trans-abdominally. The 
user moves the device around on the maternal abdomen and, using 2D image processing to 
locate the amniotic fluid areas, the device gives feedback to the user about where to acquire 
15 the 3D image data sets. The user acquires one or more 3D image data sets covering all of the 
amniotic fluid in the uterus and the data sets are then stored in the device or transferred to a 
host computer. 

The 3D ultrasound device is configured to acquire the 3D image data sets in two 
formats. The first format is a collection of two-dimensional scanplanes, each scanplane 
20 being separated from the other and representing a portion of the uterus being scanned. Each 
scanplane is formed from one-dimensional ultrasound A-lines confined within the limits of 
the 2D scanplane. The 3D data sets is then represented as a 3D array of 2D scanplanes. The 
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3D array of 2D scanplanes is an assembly of scanplanes, and may be assembled into a 
translational array, a wedge array, or a rotatational array. 

Alternatively, the 3D ultrasound device is configured to acquire the 3D image data 
sets from one-dimensional ultrasound A-lines distributed in 3D space of the uterus to form a 
5 3D scancone of 3D-distributed scanline. The 3D scancone is not an assembly of 2D 
scanplanes. 

The 3D image datasets, either as discrete scanplanes or 3D distributed scanlines, are 
then subjected to image enhancement and analysis processes. The processes are either 
implemented on the device itself or is implemented on the host computer. Alternatively, the 
10 processes can also be implemented on a server or other computer to which the 3D ultrasound 
data sets are transferred. 

In a preferred image enhancement process, each 2D image in the 3D dataset is first 
enhanced using non-linear filters by an image pre-filtering step. The image pre-filtering step 
includes an image-smoothing step to reduce image noise followed by an image-sharpening 
1 5 step to obtain maximum contrast between organ wall boundaries. 

A second process includes subjecting the resulting image of the first process to a 
location method to identify initial edge points between amniotic fluid and other fetal or 
maternal structures. The location method automatically determines the leading and trailing 
regions of wall locations along an A-mode one-dimensional scan line. 

20 A third process includes subjecting the image of the first process to an intensity-based 

segmentation process where dark pixels (representing fluid) are automatically separated from 
bright pixels (representing tissue and other structures). 
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In a fourth process, the images resulting from the second and third step are combined 
to result in a single image representing likely amniotic fluid regions. 

In a fifth process, the combined image is cleaned to make the output image smooth 
and to remove extraneous structures such as the fetal head and the fetal bladder. 
5 In a sixth process, boundary line contours are placed on each 2D image. Thereafter, 

the method then calculates the total 3D volume of amniotic fluid in the uterus. 

In cases in which uteruses are too large to fit in a single 3D array of 2D scanplanes or 
a single 3D scancone of 3D distributed scanlines, especially as occurs during the second and 
third trimester of pregnancy, preferred alternate embodiments of the invention allow for 
10 acquiring at least two 3D data sets, preferably four, each 3D data set having at least a partial 
ultrasonic view of the uterus, each partial view obtained from a different anatomical site of 
the patient. 

In one embodiment a 3D array of 2D scanplanes is assembled such that the 3D array 
presents a composite image of the uterus that displays the amniotic fluid regions to provide 

15 the basis for calculation of amniotic fluid volumes. In a preferred alternate embodiment, the 
user acquires the 3D data sets in quarter sections of the uterus when the patient is in a supine 
position. In this 4-quadrant supine procedure, four image cones of data are acquired near the 
midpoint of each uterine quadrant at substantially equally spaced intervals between quadrant 
centers. Image processing as outlined above is conducted for each quadrant image, 

20 segmenting on the darker pixels or voxels associated with amniotic fluid. Correcting 
algorithms are applied to compensate for any quadrant-to-quadrant image cone overlap by 
registering and fixing one quadrant's image to another. The result is a fixed 3D mosaic 
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image of the uterus and the amniotic fluid volumes or regions in the uterus from the four 
separate image cones. 

Similarly, in another preferred alternate embodiment, the user acquires one or more 
3D image data sets of quarter sections of the uterus when the patient is in a lateral position. 
5 In this multi-image cone lateral procedure, each image cones of data are acquired along a 
lateral line of substantially equally spaced intervals. Each image cone are subjected to the 
image processing as outlined above, with emphasis given to segmenting on the darker pixels 
or voxels associated with amniotic fluid. Scanplanes showing common pixel or voxel 
overlaps are registered into a common coordinate system along the lateral line. Correcting 
10 algorithms are applied to compensate for any image cone overlap along the lateral line. The 
result is a fixed 3D mosaic image of the uterus and the amniotic fluid volumes or regions in 
the uterus from the four separate image cones. 

In yet other preferred embodiments, at least two 3D scancone of 3D distributed 
scanlines are acquired at different anatomical sites, image processed, registered and fused 
15 into a 3D mosaic image composite. Amniotic fluid volumes are then calculated. 

The system and method further provides an automatic method to detect and correct 
for any contribution the fetal head provides to the amniotic fluid volume. 

Brief Description of the Drawings 
FIGURE 1 is a side view of a microprocessor-controlled, hand-held ultrasound 

20 transceiver; 

FIGURE 2A is a is depiction of the hand-held transceiver in use for scanning a patient; 
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FIGURE 2B is a perspective view of the hand-held transceiver device sitting in a 
communication cradle; 

FIGURE 3 depicts a schematic view of a plurality of transceivers in connection with a 

server; 

FIGURE 4 depicts a schematic view of a plurality of transceivers in connection with a 
server over a network; 

FIGURE 5 A a graphical representation of a plurality of scan lines forming a single scan 

plane; 

FIGURE 5B is a graphical representation of a plurality of scanplanes forming a three- 
dimensional array having a substantially conic shape; 

FIGURE 5C is a graphical representation of a plurality of 3D distributed scanlines 
emanating from the transceiver forming a scancone; 

FIGURE 6 is a depiction of the hand-held transceiver placed laterally on a patient trans- 
abdominally to transmit ultrasound and receive ultrasound echoes for processing to determine 
amniotic fluid volumes; 

FIGURE 7 shows a block diagram overview of the two-dimensional and three- 
dimensional Input, Image Enhancement, Intensity-Based Segmentation, Edge-Based 
Segmentation, Combine, Polish , Output, and Compute algorithms to visualize and determine 
the volume or area of amniotic fluid; 

FIGURE 8 A depicts the sub-algorithms of Image Enhancement; 

FIGURE 8B depicts the sub-algorithms of Intensity-Based Segmentation; 
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FIGURE 8C depicts the sub-algorithms of Edge-Based Segmentation; 

FIGURE 8D depicts the sub-algorithms of the Polish algorithm, including Close, Open, 
Remove Deep Regions, and Remove Fetal Head Regions; 

FIGURE 8E depicts the sub-algorithms of the Remove Fetal Head Regions sub- 
5 algorithm; 

FIGURE 8F depicts the sub-algorithms of the Hough Transform sub-algorithm; 

FIGURE 9 depicts the operation of a circular Hough transform algorithm; 

FIGURE 10 shows results of sequentially applying the algorithm steps on a sample 

image; 

10 FIGURE 1 1 illustrates a set of intermediate images of the fetal head detection process; 

FIGURE 12 presents a 4-panel series of sonographer amniotic fluid pocket outlines and 
the algorithm output amniotic fluid pocket outlines; 

FIGURE 13 illustrates a 4-quadrant supine procedure to acquire multiple image cones; 

FIGURE 14 illustrates an in-line lateral line procedure to acquire multiple image cones; 

15 FIGURE 15 is a block diagram overview of the rigid registration and correcting 

algorithms used in processing multiple image cone data sets; 

FIGURE 16 is a block diagram of the steps in the rigid registration algorithm; 

FIGURE 17A is an example image showing a first view of a fixed scanplane; 

FIGURE 17B is an example image showing a second view view of a moving scanplane 
20 having some voxels in common with the first scanplane; 

FIGURE 17C is a composite image of the first (fixed) and second (moving) images; 
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FIGURE 18A is an example image showing a first view of a fixed scanplane; 

FIGURE 18B is an example image showing a second view of a moving scanplane 
having some voxels in common with the first view and a third view; 

FIGURE 18C is a third view of a moving scanplane having some voxels in common 
5 with the second view; and 

FIGURE 18D is a composite image of the first (fixed), second (moving), and third 
(moving) views. 

Detailed Description of the Preferred Embodiment 

The preferred portable embodiment of the ultrasound transceiver of the amniotic fluid 
10 volume measuring system are shown in FIGURES 1-4. The transceiver 10 includes a handle 
12 having a trigger 14 and a top button 16, a transceiver housing 18 attached to the handle 12, 
and a transceiver dome 20. A display 24 for user interaction is attached to the transceiver 
housing 18 at an end opposite the transceiver dome 20. Housed within the transceiver 10 is a 
single element transducer (not shown) that converts ultrasound waves to electrical signals. The 
15 transceiver 10 is held in position against the body of a patient by a user for image acquisition 
and signal processing. In operation, the transceiver 10 transmits a radio frequency ultrasound 
signal at substantially 3.7 MHz to the body and then receives a returning echo signal. To 
accommodate different patients having a variable range of obesity, the transceiver 10 can be 
adjusted to transmit a range of probing ultrasound energy from approximately 2 MHz to 
20 approximately 10 MHz radio frequencies. 

The top button 16 selects for different acquisition volumes. The transceiver is 
controlled by a microprocessor and software associated with the microprocessor and a digital 
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signal processor of a computer system. As used in this invention, the term "computer system" 
broadly comprises any microprocessor-based or other computer system capable of executing 
operating instructions and manipulating data, and is not limited to a traditional desktop or 
notebook computer. The display 24 presents alphanumeric or graphic data indicating the 
5 proper or optimal positioning of the transceiver 10 for initiating a series of scans. The 
transceiver 10 is configured to initiate the series of scans to obtain and present 3D images as 
either a 3D array of 2D scanplanes or as a single 3D scancone of 3D distributed scanlines. A 
suitable transceiver is the DCD372 made by Diagnostic Ultrasound. In alternate embodiments, 
the two- or three-dimensional image of a scan plane may be presented in the display 24. 

10 Although the preferred ultrasound transceiver is described above, other transceivers 

may also be used. For example, the transceiver need not be battery-operated or otherwise 
portable, need not have a top-mounted display 24, and may include many other features or 
differences. The display 24 may be a liquid crystal display (LCD), a light emitting diode 
(LED), a cathode ray tube (CRT), or any suitable display capable of presenting alphanumeric 

1 5 data or graphic images. 

FIGURE 2 A is a photograph of the hand-held transceiver 10 for scanning a patient. 
The transceiver 10 is then positioned over the patient's abdomen by a user holding the handle 
12 to place the transceiver housing 18 against the patient's abdomen. The top button 16 is 
centrally located on the handle 12. Once optimally positioned over the abdomen for scanning, 
20 the transceiver 10 transmits an ultrasound signal at substantially 3.7 MHz into the uterus. The 
transceiver 10 receives a return ultrasound echo signal emanating from the uterus and presents 
it on the display 24. 
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FIGURE 2B is a perspective view of the hand-held transceiver device sitting in a 
communication cradle. The transceiver 10 sits in a communication cradle 42 via the handle 12. 
This cradle can be connected to a standard USB port of any personal computer, enabling all the 
data on the device to be transferred to the computer and enabling new programs to be 
5 transferred into the device from the computer. 

FIGURE 3 depicts a schematic view of a plurality of transceivers in connection with a 
server. FIGURE 3, by example, depicts each transceiver 10 being used to send probing 
ultrasound radiation to a uterus of a patient and to subsequently retrieve ultrasound echoes 
returning from the uterus, convert the ultrasound echoes into digital echo signals, store the 

10 digital echo signals, and process the digital echo signals by algorithms of the invention. A user 
holds the transceiver 10 by the handle 12 to send probing ultrasound signals and to receive 
incoming ultrasound echoes. The transceiver 10 is placed in the communication cradle 42 that 
is in signal communication with a computer 52, and operates as an amniotic fluid volume 
measuring system. Two amniotic fluid volume-measuring systems are depicted as 

15 representative though fewer or more systems may be used. As used in this invention, a 
"server" can be any computer software or hardware that responds to requests or issues 
commands to or from a client. Likewise, the server may be accessible by one or more client 
computers via the Internet, or may be in communication over a LAN or other network. 

Each amniotic fluid volume measuring systems includes the transceiver 10 for 
20 acquiring data from a patient. The transceiver 10 is placed in the cradle 52 to establish signal 
communication with the computer 52. Signal communication as illustrated is by a wired 
connection from the cradle 42 to the computer 52. Signal communication between the 
transceiver 10 and the computer 52 may also be by wireless means, for example, infrared 
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signals or radio frequency signals. The wireless means of signal communication may occur 
between the cradle 42 and the computer 52, the transceiver 10 and the computer 52, or the 
transceiver 10 and the cradle 42. 

A preferred first embodiment of the amniotic fluid volume measuring system includes 
5 each transceiver 10 being separately used on a patient and sending signals proportionate to the 
received and acquired ultrasound echoes to the computer 52 for storage. Residing in each 
computer 52 are imaging programs having instructions to prepare and analyze a plurality of 
one dimensional (ID) images from the stored signals and transforms the plurality of ID images 
into the plurality of 2D scanplanes. The imaging programs also present 3D renderings from the 
10 plurality of 2D scanplanes. Also residing in each computer 52 are instructions to perform the 
additional ultrasound image enhancement procedures, including instructions to implement the 
image processing algorithms. 

A preferred second embodiment of the amniotic fluid volume measuring system is 
similar to the first embodiment, but the imaging programs and the instructions to perform the 

15 additional ultrasound enhancement procedures are located on the server 56. Each computer 52 
from each amniotic fluid volume measuring system receives the acquired signals from the 
transceiver 10 via the cradle 51 and stores the signals in the memory of the computer 52. The 
computer 52 subsequently retrieves the imaging programs and the instructions to perform the 
additional ultrasound enhancement procedures from the server 56. Thereafter, each computer 

20 52 prepares the ID images, 2D images, 3D renderings, and enhanced images from the retrieved 
imaging and ultrasound enhancement procedures. Results from the data analysis procedures 
are sent to the server 56 for storage. 



- 15 



Black Lowe & Graham f 



25315 

PATENT TRADEMARK OFFICE 



DXUC-1-1020AP 



816 Second Avenue 
Seattle, Washington 98104 
206.381.3300 • F: 206.381.3301 



A preferred third embodiment of the amniotic fluid volume measuring system is similar 
to the first and second embodiments, but the imaging programs and the instructions to perform 
the additional ultrasound enhancement procedures are located on the server 56 and executed on 
the server 56. Each computer 52 from each amniotic fluid volume measuring system receives 
5 the acquired signals from the transceiver 10 and via the cradle 51 sends the acquired signals in 
the memory of the computer 52. The computer 52 subsequently sends the stored signals to the 
server 56. In the server 56, the imaging programs and the instructions to perform the additional 
ultrasound enhancement procedures are executed to prepare the ID images, 2D images, 3D 
renderings, and enhanced images from the server 56 stored signals. Results from the data 
10 analysis procedures are kept on the server 56, or alternatively, sent to the computer 52. 

FIGURE 4 is a schematic view of a plurality of amniotic fluid measuring systems 
connected to a server over the Internet or other network 64. FIGURE 4 represents any of the 
first, second, or third embodiments of the invention advantageously deployed to other servers 
and computer systems through connections via the network. 

1 5 FIGURE 5 A a graphical representation of a plurality of scan lines forming a single scan 

plane. FIGURE 5A illustrates how ultrasound signals are used to make analyzable images, 
more specifically how a series of one-dimensional (ID) scanlines are used to produce a two- 
dimensional (2D) image. The ID and 2D operational aspects of the single element transducer 
housed in the transceiver 10 is seen as it rotates mechanically about an angle <(). A scanline 214 

20 of length r migrates between a first limiting position 218 and a second limiting position 222 as 
determined by the value of the angle <|>, creating a fan-like 2D scanplane 210. In one preferred 
form, the transceiver 10 operates substantially at 3.7 MHz frequency and creates an 
approximately 18 cm deep scan line 214 and migrates within the angle <j> having an angle of 
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approximately 0.027 radians. A first motor tilts the transducer approximately 60° clockwise 
and then counterclockwise forming the fan-like 2D scanplane presenting an approximate 120° 
2D sector image. A plurality of scanlines, each scanline substantially equivalent to scanline 
214 is recorded, between the first limiting position 218 and the second limiting position 222 
5 formed by the unique tilt angle ((>. The plurality of scanlines between the two extremes forms a 
scanplane 210. In the preferred embodiment, each scanplane contains 77 scan lines, although 
the number of lines can vary within the scope of this invention. The tilt angle <|) sweeps 
through angles approximately between -60° and +60° for a total arc of approximately 120°. 

FIGURE 5B is a graphical representation of a plurality of scanplanes forming a three- 
10 dimensional array (3D) 240 having a substantially conic shape. FIGURE 5B illustrates how a 
3D rendering is obtained from the plurality of 2D scanplanes. Within each scanplane 210 are 
the plurality of scanlines, each scanline equivalent to the scanline 214 and sharing a common 
rotational angle 0. In the preferred embodiment, each scanplane contains 77 scan lines, 
although the number of lines can vary within the scope of this invention. Each 2D sector 
15 image scanplane 210 with tilt angle <j) and range r (equivalent to the scanline 214) collectively 
forms a 3D conic array 240 with rotation angle 0. After gathering the 2D sector image, a 
second motor rotates the transducer between 3.75° or 7.5° to gather the next 120° sector image. 
This process is repeated until the transducer is rotated through 180°, resulting in the cone- 
shaped 3D conic array 240 data set with 24 planes rotationally assembled in the preferred 
20 embodiment. The conic array could have fewer or more planes rotationally assembled. For 
example, preferred alternate embodiments of the conic array could include at least two 
scanplanes, or a range of scanplanes from 2 to 48 scanplanes. The upper range of the 
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scanplanes can be greater than 48 scanplanes. The tilt angle § indicates the tilt of the scanline 
from the centerline in 2D sector image, and the rotation angle 0, identifies the particular 
rotation plane the sector image lies in. Therefore, any point in this 3D data set can be isolated 
using coordinates expressed as three parameters, P(r,<f>,0). 

5 As the scanlines are transmitted and received, the returning echoes are interpreted as 

analog electrical signals by a transducer, converted to digital signals by an analog-to-digital 
converter, and conveyed to the digital signal processor of the computer system for storage and 
analysis to determine the locations of the amniotic fluid walls. The computer system is 
representationally depicted in FIGURES 3 and 4 and includes a microprocessor, random access 
1 0 memory (RAM), or other memory for storing processing instructions and data generated by the 
transceiver 10. 

FIGURE 5C is a graphical representation of a plurality of 3D-distributed scanlines 
emanating from the transceiver 10 forming a scancone 300. The scancone 300 is formed by a 
plurality of 3D distributed scanlines that comprises a plurality of internal and peripheral 

15 scanlines. The scanlines are one-dimensional ultrasound A-lines that emanate from the 
tranciever 10 at different coordinate directions, that taken as an aggregate, from a conic shape. 
The 3D-distributed A-lines (scanlines) are not necessarily confined within a scanplane, but 
instead are directed to sweep throughout the internal and along the periphery of the scancone 
300. The 3D-distributed scanlines not only would occupy a given scanplane in a 3D array of 

20 2D scanplanes, but also the inter-scanplane spaces, from the conic axis to and including the 
conic periphery. The transceiver 10 shows the same illustrated features from FIGURE 1 , but is 
configured to distribute the ultrasound A-lines throughout 3D space in different coordinate 
directions to form the scancone 300. 
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The internal scanlines are represented by scanlines 312A-C. The number and location 
of the internal scanlines emanating from the transceiver 10 is the number of internal scanlines 
needed to be distributed within the scancone 300, at different positional coordinates, to 
sufficiently visualize structures or images within the scancone 300. The internal scanlines are 
5 not peripheral scanlines. The peripheral scanlines are represented by scanlines 314A-F and 
occupy the conic periphery, thus representing the peripheral limits of the scancone 300. 

FIGURE 6 is a depiction of the hand-held transceiver placed on a patient trans- 
abdominal^ to transmit probing ultrasound and receive ultrasound echoes for processing to 
determine amniotic fluid volumes. The transceiver 10 is held by the handle 12 to position 
10 over a patient to measure the volume of amniotic fluid in an amniotic sac over a baby. A 
plurality of axes for describing the orientation of the baby, the amniotic sac, and mother is 
illustrated. The plurality of axes includes a vertical axis depicted on the line L(R) - L(L) for 
left and right orientations, a horizontal axis LI - LS for inferior and superior orientations, and 
a depth axis LA - LP for anterior and posterior orientations. 

1 5 FIGURE 6 is representative of a preferred data acquisition protocol used for amniotic 

fluid volume determination. In this protocol, the transceiver 10 is the hand-held 3D 
ultrasound device (for example, model DCD372 from Diagnostic Ultrasound) and is used to 
image the uterus trans-abdominally. Initially during the targeting phase, the patient is in a 
supine position and the device is operated in a 2D continuous acquisition mode. A 2D 

20 continuous mode is where the data is continuously acquired in 2D and presented as a 
scanplane similar to the scanplane 210 on the display 24 while an operator physically moves 
the transceiver 10. An operator moves the transceiver 10 around on the maternal abdomen 
and the presses the trigger 14 of the transceiver 10 and continuously acquires real-time 
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feedback presented in 2D on the display 24. Amniotic fluid, where present, visually appears 
as dark regions along with an alphanumeric indication of amniotic fluid area (for example, in 
cm 2 ) on the display 24. Based on this real-time information in terms of the relative position 
of the transceiver 10 to the fetus, the operator decides which side of the uterus has more 
5 amniotic fluid by the presentation on the display 24. The side having more amniotic fluid 
presents as regions having larger darker regions on the display 24. Accordingly, the side 
displaying a large dark region registers greater alphanumeric area while the side with less 
fluid shows displays smaller dark regions and proportionately registers smaller alphanumeric 
area on the display 24. While amniotic fluid is present throughout the uterus, its distribution 
10 in the uterus depends upon where and how the fetus is positioned within the uterus. There is 
usually less amniotic fluid around the fetus's spine and back and more amniotic fluid in front 
of its abdomen and around the limbs. 

Based on fetal position information acquired from data gathered under continuous 
acquisition mode, the patient is placed in a lateral recumbent position such that the fetus is 
15 displaced towards the ground creating a large pocket of amniotic fluid close to abdominal 
surface where the transceiver 10 can be placed as shown in FIGURE 6. For example, if large 
fluid pockets are found on the right side of the patient, the patient is asked to turn with the 
left side down and if large fluid pockets are found on the left side, the patient is asked to turn 
with the right side down. 

20 After the patient has been placed in the desired position, the transceiver 10 is again 

operated in the 2D continuous acquisition mode and is moved around on the lateral surface of 
the patient's abdomen. The operator finds the location that shows the largest amniotic fluid 
area based on acquiring the largest dark region imaged and the largest alphanumeric value 
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displayed on the display 24. At the lateral abdominal location providing the largest dark 
region, the transceiver 10 is held in a fixed position, the trigger 14 is released to acquire a 3D 
image comprising a set of arrayed scanplanes. The 3D image presents a rotational array of 
the scanplanes 210 similar to the 3D array 240. 

5 In a preferred alternate data acquisition protocol, the operator can reposition the 

transceiver 10 to different abdominal locations to acquire new 3D images comprised of 
different scanplane arrays similar to the 3D array 240. Multiple scan cones obtained from 
different positions provide the operator the ability to image the entire amniotic fluid region 
from different view points. In the case of a single image cone being too small to 
10 accommodate a large AFV measurement, obtaining multiple 3D array 240 image cones 
ensures that the total volume of large AFV regions is determined. Multiple 3D images may 
also be acquired by pressing the top bottom 16 to select multiple conic arrays similar to the 
3D array 240. 

Depending on the position of the fetus relative to the location of the transceiver 10, a 
15 single image scan may present an underestimated volume of AFV due to amniotic fluid 

pockets that remain hidden behind the limbs of the fetus. The hidden amniotic fluid pockets 

present as unquantifiable shadow-regions. 

To guard against underestimating AFV, repeated positioning the transceiver 10 and 

rescanning can be done to obtain more than one ultrasound view to maximize detection of 
20 amniotic fluid pockets. Repositioning and rescanning provides multiple views as a plurality of 

the 3D arrays 240 images cones. Acquiring multiple images cones improves the probability of 

obtaining initial estimates of AFV that otherwise could remain undetected and un-quantified in 

a single scan. 
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In an alternative scan protocol, the user determines and scans at only one location on 
the entire abdomen that shows the maximum amniotic fluid area while the patient is the 
supine position. As before, when the user presses the top button 16, 2D scanplane images 
equivalent to the scanplane 210 are continuously acquired and the amniotic fluid area on 
5 every image is automatically computed. The user selects one location that shows the 
maximum amniotic fluid area. At this location, as the user releases the scan button, a full 3D 
data cone is acquired and stored in the device's memory. 

FIGURE 7 shows a block diagram overview the image enhancement, segmentation, 
and polishing algorithms of the amniotic fluid volume measuring system. The enhancement 

10 segmentation, and polishing algorithms are applied to each scanplane 210 or to the entire scan 
cone 240 to automatically obtain amniotic fluid regions. For scanplanes substantially 
equivalent to scanplane 210, the algorithms are expressed in two-dimensional terms and use 
formulas to convert scanplane pixels (picture elements) into area units. For the scan cones 
substantially equivalent to the 3D conic array 240, the algorithms are expressed in three- 

1 5 dimensional terms and use formulas to convert voxels (volume elements) into volume units. 

The algorithms expressed in 2D terms are used during the targeting phase where the 
operator trans-abdominally positions and repositions the transceiver 10 to obtain real-time 
feedback about the amniotic fluid area in each scanplane. The algorithms expressed in 3D 
terms are used to obtain the total amniotic fluid volume computed from the voxels contained 
20 within the calculated amniotic fluid regions in the 3D conic array 240. 

FIGURE 7 represents an overview of a preferred method of the invention and includes 
a sequence of algorithms, many of which have sub-algorithms described in more specific detail 
in FIGURES 8A-F. FIGURE 7 begins with inputting data of an unprocessed image at step 
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410. After unprocessed image data 410 is entered (e.g., read from memory, scanned, or 
otherwise acquired), it is automatically subjected to an image enhancement algorithm 418 that 
reduces the noise in the data (including speckle noise) using one or more equations while 
preserving the salient edges on the image using one or more additional equations. Next, the 
5 enhanced images are segmented by two different methods whose results are eventually 
combined. A first segmentation method applies an intensity-based segmentation algorithm 422 
that determines all pixels that are potentially fluid pixels based on their intensities. A second 
segmentation method applies an edge-based segmentation algorithm 438 that relies on 
detecting the fluid and tissue interfaces. The images obtained by the first segmentation 

10 algorithm 422 and the images obtained by the second segmentation algorithm 438 are brought 
together via a combination algorithm 442 to provide a substantially segmented image. The 
segmented image obtained from the combination algorithm 442 are then subjected to a 
polishing algorithm 464 in which the segmented image is cleaned-up by filling gaps with pixels 
and removing unlikely regions. The image obtained from the polishing algorithm 464 is 

15 outputted 480 for calculation of areas and volumes of segmented regions-of-interest. Finally 
the area or the volume of the segmented region-of-interest is computed 484 by multiplying 
pixels by a first resolution factor to obtain area, or voxels by a second resolution factor to 
obtain volume. For example, for pixels having a size of 0.8mm by 0.8 mm, the first resolution 
or conversion factor for pixel area is equivalent to 0.64 mm 2 , and the second resolution or 

20 conversion factor for voxel volume is equivalent to 0.512 mm 3 . Different unit lengths for 
pixels and voxels may be assigned, with a proportional change in pixel area and voxel volume 
conversion factors. 
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The enhancement, segmentation and polishing algorithms depicted in FIGURE 7 for 
measuring amniotic fluid areas or volumes are not limited to scanplanes assembled into 
rotational arrays equivalent to the 3D array 240. As additional examples, the enhancement, 
segmentation and polishing algorithms depicted in FIGURE 7 apply to translation arrays and 
wedge arrays. Translation arrays are substantially rectilinear image plane slices from 
incrementally repositioned ultrasound transceivers that are configured to acquire ultrasound 
rectilinear scanplanes separated by regular or irregular rectilinear spaces. The translation 
arrays can be made from transceivers configured to advance incrementally, or may be hand- 
positioned incrementally by an operator. The operator obtains a wedge array from ultrasound 
transceivers configured to acquire wedge-shaped scanplanes separated by regular or irregular 
angular spaces, and either mechanistically advanced or hand-tilted incrementally. Any number 
of scanplanes can be either translationally assembled or wedge-assembled ranges, but 
preferably in ranges greater than 2 scanplanes. 

Other preferred embodiments of the enhancement, segmentation and polishing 
algorithms depicted in FIGURE 7 may be applied to images formed by line arrays, either spiral 
distributed or reconstructed random-lines. The line arrays are defined using points identified 
by the coordinates expressed by the three parameters, P(r, </>,&), where the values or r , <f> , and 
0 can vary. 

The enhancement, segmentation and polishing algorithms depicted in FIGURE 7 are 
not limited to ultrasound applications but may be employed in other imaging technologies 
utilizing scanplane arrays or individual scanplanes. For example, biological-based and non- 
biological-based images acquired using infrared, visible light, ultraviolet light, microwave, x- 
ray computed tomography, magnetic resonance, gamma rays, and positron emission are images 
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suitable for the algorithms depicted in FIGURE 7. Furthermore, the algorithms depicted in 
FIGURE 7 can be applied to facsimile transmitted images and documents. 

FIGURES 8A-E depict expanded details of the preferred embodiments of 
enhancement, segmentation, and polishing algorithms described in Figure 7. Each of the 
5 following greater detailed algorithms are either implemented on the transceiver 10 itself or are 
implemented on the host computer 52 or on the server 56 computer to which the ultrasound 
data is transferred. 

FIGURE 8A depicts the sub-algorithms of Image Enhancement. The sub-algorithms 
include a heat filter 514 to reduce noise and a shock niter 518 to sharpen edges. A 

10 combination of the heat and shock filters works very well at reducing noise and sharpening the 
data while preserving the significant discontinuities. First, the noisy signal is filtered using a 
ID heat filter (Equation El below), which results in the reduction of noise and smoothing of 
edges. This step is followed by a shock-filtering step 518 (Equation E2 below), which results in 
the sharpening of the blurred signal. Noise reduction and edge sharpening is achieved by 

15 application of the following equations E1-E2. The algorithm of the heat filter 514 uses a heat 
equation El. The heat equation El in partial differential equation (PDE) form for image 
processing is expressed as: 

du_6^u tfu 

where u is the image being processed. The image u is 2D, and is comprised of an array 
20 of pixels arranged in rows along the x-axis, and an array of pixels arranged in columns along 
the y-axis. The pixel intensity of each pixel in the image u has an initial input image pixel 
intensity (I) defined as u 0 = I . The value of / depends on the application, and commonly 
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occurs within ranges consistent with the application. For example, / can be as low as 0 to 1, or 
occupy middle ranges between 0 to 127 or 0 to 512. Similarly, / may have values occupying 
higher ranges of 0 to 1024 and 0 to 4096, or greater. The heat equation El results in a 
smoothing of the image and is equivalent to the Gaussian filtering of the image. The larger the 
5 number of iterations that it is applied for the more the input image is smoothed or blurred and 
the more the noise that is reduced. 

The shock filter 518 is a PDE used to sharpen images as detailed below. The two 
dimensional shock filter E2 is expressed as: 

^-F(t(u)pu\ 9 E2 

10 where u is the image processed whose initial value is the input image pixel intensity (I): u 0 = I 
where the £(u) term is the Laplacian of the image u, F is a function of the Laplacian, and ||Vw|| 
is the 2D gradient magnitude of image intensity defined by equation E3. 

where 

15 u 2 x = the square of the partial derivative of the pixel intensity (u) along the x- 

axis, 

u 2 y = the square of the partial derivative of the pixel intensity (u) along the y- 
axis, 

the Laplacian £(u) of the image, u, is expressed in equation E4 as 

20 = U xx U x 2 + 2U xy U x U y + U yy U y E 4 
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where equation E4 relates to equation El as follows: 



du 

ux is the first partial derivative — of u along the x-axis, 

dx 



du 

u y is the first partial derivative — of u along the y-axis, 

dy 



du 

u 2 is the square of the first partial derivative — of u along the x-axis, 

dx 



du 

u is the square of the first partial derivative — of u along the y-axis, 

dy 



d 2 u 

uxx is the second partial derivative — r-of u along the x-axis, 

dx 



dy 2 



u yy is the second partial derivative of u along the y-axis, 



du 

Uxy 'is cross multiple first partial derivative of u along the x and y 

dxdy 

axes, and 

the sign of the function F modifies the Laplacian by the image gradient values 
selected to avoid placing spurious edges at points with small gradient values: 



F(£(u)) = 1, if £(u) > 0 and[|Vw|| > t 

= -l,if £(u)<0 and||Vw||>/ 
= 0, otherwise 

where / is a threshold on the pixel gradient value ||Vw| 
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The combination of heat filtering and shock filtering produces an enhanced image 
ready to undergo the intensity-based and edge-based segmentation algorithms as discussed 
below. 

FIGURE 8B depicts the sub-algorithms of Intensity-Based Segmentation (step 422 in 
5 FIGURE 7). The intensity-based segmentation step 422 uses a "k-means" intensity 
clustering 522 technique where the enhanced image is subjected to a categorizing "k-means" 
clustering algorithm. The "k-means" algorithm categorizes pixel intensities into white, gray, 
and black pixel groups. Given the number of desired clusters or groups of intensities (k), the 
k-means algorithm is an iterative algorithm comprising four steps: 

10 1. Initially determine or categorize cluster boundaries by defining a 

minimum and a maximum pixel intensity value for every white, gray, or black pixels into 
groups or k-clusters that are equally spaced in the entire intensity range. 

2. Assign each pixel to one of the white, gray or black k-clusters based on 

the currently set cluster boundaries. 

15 3. Calculate a mean intensity for each pixel intensity k-cluster or group based 

on the current assignment of pixels into the different k-clusters. The calculated mean 
intensity is defined as a cluster center. Thereafter, new cluster boundaries are determined as 
mid points between cluster centers. 

4. Determine if the cluster boundaries significantly change locations from 

20 their previous values. Should the cluster boundaries change significantly from their previous 
values, iterate back to step 2, until the cluster centers do not change significantly between 
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iterations. Visually, the clustering process is manifest by the segmented image and repeated 
iterations continue until the segmented image does not change between the iterations. 

The pixels in the cluster having the lowest intensity value - the darkest cluster - are 
defined as pixels associated with amniotic fluid. For the 2D algorithm, each image is 
clustered independently of the neighboring images. For the 3D algorithm, the entire volume 
is clustered together. To make this step faster, pixels are sampled at 2 or any multiple 
sampling rate factors before determining the cluster boundaries. The cluster boundaries 
determined from the down-sampled data are then applied to the entire data. 

FIGURE 8C depicts the sub-aigorithms of Edge-Based Segmentation (step 438 in 
FIGURE 7) and uses a sequence of four sub-algorithms. The sequence includes a spatial 
gradients 526 algorithm, a hysteresis threshold 530 algorithm, a Region-of-Interest (ROI) 
534 algorithm, and a matching edges filter 538 algorithm. 

The spatial gradient 526 computes the x-directional and y-directional spatial gradients 
of the enhanced image. The Hysteresis threshold 530 algorithm detects salient edges. Once 
the edges are detected, the regions defined by the edges are selected by a user employing the 
ROI 534 algorithm to select regions-of-interest deemed relevant for analysis. 

Since the enhanced image has very sharp transitions, the edge points can be easily 
determined by taking x- and y- derivatives using backward differences along x- and y- 
directions. The pixel gradient magnitude ||W|| is then computed from the x- and y-derivative 
image in equation E5 as: 




E5 



Where / x = the square of x-derivative of intensity; and 
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fy = the square of y-derivative of intensity along the y-axis. 

Significant edge points are then determined by thresholding the gradient magnitudes 
using a hysteresis thresholding operation. Other thresholding methods could also be used. In 
hysteresis thresholding 530, two threshold values, a lower threshold and a higher threshold, 
5 are used. First, the image is thresholded at the lower threshold value and a connected 
component labeling is carried out on the resulting image. Next, each connected edge 
component is preserved which has at least one edge pixel having a gradient magnitude 
greater than the upper threshold. This kind of thresholding scheme is good at retaining long 
connected edges that have one or more high gradient points. 

10 In the preferred embodiment, the two thresholds are automatically estimated. The 

upper gradient threshold is estimated at a value such that at most 97% of the image pixels are 
marked as non-edges. The lower threshold is set at 50% of the value of the upper threshold. 
These percentages could be different in different implementations. Next, edge points that lie 
within a desired region-of-interest are selected 534. This region of interest selection 534 

15 excludes points lying at the image boundaries and points lying too close to or too far from the 
transceiver 10. Finally, the matching edge filter 538 is applied to remove outlier edge points 
and fill in the area between the matching edge points. 

The edge-matching algorithm 538 is applied to establish valid boundary edges and 
remove spurious edges while filling the regions between boundary edges. Edge points on an 
20 image have a directional component indicating the direction of the gradient. Pixels in 
scanlines crossing a boundary edge location will exhibit two gradient transitions depending 
on the pixel intensity directionality. Each gradient transition is given a positive or negative 
value depending on the pixel intensity directionality. For example, if the scanline approaches 
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an echo reflective bright wall from a darker region, then an ascending transition is 
established as the pixel intensity gradient increases to a maximum value, i.e., as the transition 
ascends from a dark region to a bright region. The ascending transition is given a positive 
numerical value. Similarly, as the scanline recedes from the echo reflective wall, a 
5 descending transition is established as the pixel intensity gradient decreases to or approaches 
a minimum value. The descending transition is given a negative numerical value. 

Valid boundary edges are those that exhibit ascending and descending pixel intensity 
gradients, or equivalently, exhibit paired or matched positive and negative numerical values. 
The valid boundary edges are retained in the image. Spurious or invalid boundary edges do 
10 not exhibit paired ascending-descending pixel intensity gradients, i.e., do not exhibit paired 
or matched positive and negative numerical values. The spurious boundary edges are 
removed from the image. 

For amniotic fluid volume related applications, most edge points for amniotic fluid 
surround a dark, closed region, with directions pointing inwards towards the center of the 

15 region. Thus, for a convex-shaped region, the direction of a gradient for any edge point, the 
edge point having a gradient direction approximately opposite to the current point represents 
the matching edge point. Those edge points exhibiting an assigned positive and negative 
value are kept as valid edge points on the image because the negative value is paired with its 
positive value counterpart. Similarly, those edge point candidates having unmatched values, 

20 i.e., those edge point candidates not having a negative-positive value pair, are deemed not to 
be true or valid edge points and are discarded from the image. 

The matching edge point algorithm 538 delineates edge points not lying on the 
boundary for removal from the desired dark regions. Thereafter, the region between any two 
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matching edge points is filled in with non-zero pixels to establish edge-based segmentation. 
In a preferred embodiment of the invention, only edge points whose directions are primarily 
oriented co-linearly with the scanline are sought to permit the detection of matching front 
wall and back wall pairs. 

5 Returning to FIGURE 7, once Intensity-Based 422 and Edge-Based Segmentation 

438 is completed, both segmentation methods use a combining step that combines the results 
of intensity-based segmentation 422 step and the edge-based segmentation 438 step using an 
AND Operator of Images 442. The AND Operator of Images 442 is achieved by a pixel- 
wise Boolean AND operator 442 step to produce a segmented image by computing the pixel 

10 intersection of two images. The Boolean AND operation 442 represents the pixels as binary 
numbers and the corresponding assignment of an assigned intersection value as a binary 
number 1 or 0 by the combination of any two pixels. For example, consider any two pixels, 
say pixel A and pixels, which can have a 1 or 0 as assigned values. If pixel A 's value is 1, and 
pixel B 's value is 1, the assigned intersection value of pixel A and pixel B is 1. If the binary 

15 value of pixel A and pixel B are both 0, or if either pixel A or pixel B is 0, then the assigned 
intersection value of pixel A and pixel B is 0. The Boolean AND operation 542 takes the 
binary any two digital images as input, and outputs a third image with the pixel values made 
equivalent to the intersection of the two input images. 

Upon completion of the AND Operator of Images 442 algorithm, the polish 464 
20 algorithm of FIGURE 7 is comprised of multiple sub-algorithms. FIGURE 8D depicts the 
sub-algorithms of the Polish 464 algorithm, including a Close 546 algorithm, an Open 550 
algorithm, a Remove Deep Regions 554 algorithm, and a Remove Fetal Head Regions 560 
algorithm. 
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Closing and opening algorithms are operations that process images based on the 
knowledge of the shape of objects contained on a black and white image, where white 
represents foreground regions and black represents background regions. Closing serves to 
remove background features on the image that are smaller than a specified size. Opening 
5 serves to remove foreground features on the image that are smaller than a specified size. The 
size of the features to be removed is specified as an input to these operations. The opening 
algorithm 550 removes unlikely amniotic fluid regions from the segmented image based on 
a-priori knowledge of the size and location of amniotic fluid pockets. 

Referring to FIGURE 8D 5 the closing 546 algorithm obtains the Apparent Amniotic 
10 Fluid Area (AAFA) or Volume (AAFV) values. The AAFA and AAFV values are 
"Apparent" and maximal because these values may contain region areas or region volumes of 
non-amniotic origin unknowingly contributing to and obscuring what otherwise would be the 
true amniotic fluid volume. For example, the AAFA and AAFV values contain the true 
amniotic volumes, and possibly as well areas or volumes due to deep tissues and undetected 
15 fetal head volumes. Thus the apparent area and volume values require correction or 
adjustments due to unknown contributions of deep tissue and of the fetal head in order to 
determine an Adjusted Amniotic Fluid Area (AdAFA) value or Volume (AdAVA) value 568. 

The AdAFA and AdAVA values obtained by the Close 546 algorithm are reduced by 
the morphological opening algorithm 550. Thereafter, the AdAFA and AdAVA values are 
20 further reduced by removing areas and volumes attributable to deep regions by using the 
Remove Deep Regions 554 algorithm. Thereafter, the polishing algorithm 464 continues by 
applying a fetal head region detection algorithm 560. 
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FIGURE 8E depicts the sub-algorithms of the Remove Fetal Head Regions sub- 
algorithm 560. The basic idea of the sub-algorithms of the fetal head detection algorithm 560 
is that the edge points that potentially represent a fetal skull are detected. Thereafter, a circle 
finding algorithm to determine the best-fitting circle to these fetal skull edges is 
5 implemented. The radii of the circles that are searched are known a priori based on the fetus' 
gestational age. The best fitting circle whose fitting metric lies above a certain pre-specified 
threshold is marked as the fetal head and the region inside this circle is the fetal head region. 
The algorithms include a gestational Age 726 input, a determine head diameter factor 730 
algorithm, a Head Edge Detection algorithm, 734, and a Hough transform procedure 736= 

10 Fetal brain tissue has substantially similar ultrasound echo qualities as presented by 

amniotic fluid. If not detected and subtracted from amniotic fluid volumes, fetal brain tissue 
volumes will be measured as part of the total amniotic fluid volumes and lead to an 
overestimation and false diagnosis of oligo or poly-hyraminotic conditions. Thus detecting 
fetal head position, measuring fetal brain matter volumes, and deducting the fetal brain 

15 matter volumes from the amniotic fluid volumes to obtain a corrected amniotic fluid volume 
serves to establish accurately measure amniotic fluid volumes. 

The gestational age input 726 begins the fetal head detection algorithm 560 and uses a 
head dimension table to obtain ranges of head bi-parietal diameters (BPD) to search for (e.g., 
30 week gestational age corresponds to a 6 cm head diameter). The head diameter range is 
20 input to both the Head Edge Detection, 734, and the Hough Transform, 736., The head edge 
detection 734 algorithm seeks out the distinctively bright ultrasound echoes from the anterior 
and posterior walls of the fetal skull while the Hough Transform algorithm, 736, finds the 
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fetal head using circular shapes as models for the fetal head in the Cartesian image (pre-scan 
conversion to polar form). 

Scanplanes processed by steps 522, 538, 530, are input to the head edge detection 
step 734. Applied as the first step in the fetal head detection algorithm 734 is the detection 
5 of the potential head edges from among the edges found by the matching edge filter. The 
matching edge 538 filter outputs pairs of edge points potentially belonging to front walls or 
back walls. Not all of these walls correspond to fetal head locations. The edge points 
representing the fetal head are determined using the following heuristics: 

(1) Looking along a one dimensional A-mode scan line, fetal head locations present a 
10 corresponding matching gradient in the opposing direction within a short distance 

approximately the same size as the thickness of the fetal skull. This distance is 
currently set to a value 1 cm. 

(2) The front wall and the back wall locations of the fetal head are within a range of 

diameters corresponding to the expected diameter 730 for the gestational age 726 of 
1 5 the fetus. Walls that are too close or too far are not likely to be head locations. 

(3) A majority of the pixels between the front and back wall locations of the fetal head lie 

within the minimum intensity cluster as defined by the output of the clustering 
algorithm 422. The percentage of pixels that need to be dark is currently defined to 
be 80%. 

20 The pixels found satisfying these features are then vertically dilated to 

produce a set of thick fetal head edges as the output of Head Edge Detection, 734. 

FIGURE 8F depicts the sub-algorithms of the Hough transform procedure 736. The 
sub-algorithms include a Polar Hough Transform 738 algorithm, a find maximum Hough 
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value 742 algorithm 742, and a fill circle region 746. The Polar Hough Transform algorithm 
looks for fetal head structures in polar coordinate terms by converting from Cartesian 
coordinates using a plurality of equations. The fetal head, which appears like a circle in a 3D 
scan-converted Cartesian coordinate image, has a different shape in the pre-scan converted 
polar space. The fetal head shape is expressed in terms of polar coordinate terms explained as 
follows: 

The coordinates of a circle in the Cartesian space (x,y) with center (x Q9 y Q ) and radius 
R are defined for an angle 0 are derived and defined in equation E5 as: 

x = .Rcos0 + jc o 
y = Rsin0 + y o 



In polar space, the coordinates (r,^) , with respect to the center (r 0 ,# 0 ) , are derived 
and defined in equation E6 as: 

rsin^ = RcosO + r 0 sin^ 0 
rcos^ = RsinO + r 0 cos^ 0 

=> (rsin^-r o sin^ 0 ) 2 +(rcos^-r 0 cos$ ) ) 2 =R 2 E6 



The Hough transform 736 algorithm using equations E5 and E6 attempts to find the best-fit 
circle to the edges of an image. A circle in the polar space is defined by a set of three 
parameters, (r 0 ,(f> 0 ,R) representing the center and the radius of the circle. 
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The basic idea for the Hough transform 736 is as follows. Suppose a circle is sought 
having a fixed radius (say, Rl) for which the best center of the circle is similarly sought. 
Now, every edge point on the input image lies on a potential circle whose center lays Rl 
pixels away from it. The set of potential centers themselves form a circle of radius Rl around 
5 each edge pixel. Now, drawing potential circles of radius Rl around each edge pixel, the 
point at which most circles intersect, a center of the circle that represents a best-fit circle to 
the given edge points is obtained. Therefore, each pixel in the Hough transform output 
contains a likelihood value that is simply the count of the number of circles passing through 
that point. 

10 FIGURE 9 illustrates the Hough Transform 736 algorithm for a plurality of circles with 

a fixed radius in a Cartesian coordinate system. A portion of the plurality of circles is 
represented by a first circle 804a, a second circle 804b, and a third circle 804c. A plurality of 
edge pixels are represented as gray squares and an edge pixel 808 is shown. A circle is drawn 
around each edge pixel to distinguish a center location 812 of a best-fit circle 816 passing 

15 through each edge pixel point; the point of the center location through which most such circles 
pass (shown by a gray star 812) is the center of the best-fit circle 816 presented as a thick dark 
line. The circumference of the best fit circle 816 passes substantially through is central portion 
of each edge pixel, represented as a series of squares substantially equivalent to the edge pixel 
808. 

20 This search for best fitting circles can be easily extended to circles with varying radii 

by adding one more degree of freedom - however, a discrete set of radii around the mean 
radii for a given gestational age makes the search significantly faster, as it is not necessary to 
search all possible radii. 
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The next step in the head detection algorithm is selecting or rejecting best-fit circles 
based on its likelihood, in the find maximum Hough Value 742 algorithm. The greater the 
number of circles passing through a given point in the Hough-space, the more likely it is to 
be the center of a best-fit circle. A 2D metric as a maximum Hough value 742 of the Hough 
5 transform 736 output is defined for every image in a dataset. The 3D metric is defined as the 
maximum of the 2D metrics for the entire 3D dataset. A fetal head is selected on an image 
depending on whether its 3D metric value exceeds a preset 3D threshold and also whether the 
2D metric exceeds a preset 2D threshold. The 3D threshold is currently set at 7 and the 2D 
threshold is currently set at 5. These thresholds have been determined by extensive training 
10 on images where the fetal head was known to be present or absent. 

Thereafter, the fetal head detection algorithm concludes with a fill circle region 746 
that incorporates pixels to the image within the detected circle. The fill circle region 746 
algorithm fills the inside of the best fitting polar circle. Accordingly, the fill circle region 746 
algorithm encloses and defines the area of the fetal brain tissue, permitting the area and 
15 volume to be calculated and deducted via algorithm 554 from the apparent amniotic fluid 
area and volume (AAFA or AAFV) to obtain a computation of the corrected amniotic fluid 
area or volume via algorithm 484. 

FIGURE 10 shows the results of sequentially applying the algorithm steps of 
FIGURES 7 and 8A-D on an unprocessed sample image 820 presented within the confines of 
20 a scanplane substantially equivalent to the scanplane 210. The results of applying the heat 
filter 514 and shock filter 518 in enhancing the unprocessed sample is shown in enhanced 
image 840. The result of intensity-based segmentation algorithms 522 is shown in image 850. 
The results of edge-based segmentation 438 algorithm using sub-algorithms 526, 530, 534 
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and 538 of the enhanced image 840 is shown in segmented image 858. The result of the 
combination 442 utilizing the Boolean AND images 442 algorithm is shown in image 862 
where white represents the amniotic fluid area. The result of applying the polishing 464 
algorithm employing algorithms 542, 546, 550, 554, 560, and 564 is shown in image 864, 
5 which depicts the amniotic fluid area overlaid on the unprocessed sample image 810. 

FIGURE 1 1 depicts a series of images showing the results of the above method to 
automatically detect, locate, and measure the area and volume of a fetal head using the 
algorithms outlined in FIGURES 7 and 8A-F. Beginning with an input image in polar 
coordinate form 920, the fetal head image is marked by distinctive bright echoes from the 
10 anterior and posterior walls of the fetal skull and a circular shape of the fetal head in the 
Cartesian image. The fetal head detection algorithm 734 operates on the polar coordinate 
data (i.e., pre-scan version, not yet converted to Cartesian coordinates). 

An example output of applying the head edge detection 734 algorithm to detect 
potential head edges is shown in image 930. Occupying the space between the anterior and 
15 posterior walls are dilated black pixels 932 (stacks or short lines of black pixels representing 
thick edges). An example of the polar Hough transform 738 for one actual data sample for a 
specific radius is shown in polar coordinate image 940. 

An example of the best-fit circle on real data polar data is shown in polar coordinate 
image 950 that has undergone the find maximum Hough value step 742. The polar 
20 coordinate image 950 is scan-converted to a Cartesian data in image 960 where the effects of 
finding maximum Hough value 742 algorithm are seen in Cartesian format. 

FIGURE 12 presents a 4-panel series of sonographer amniotic fluid pocket outlines 
compared to the algorithm's output in a scanplane equivalent to scanplane 210. The top two 
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panels depict the sonographer's outlines of amniotic fluid pockets obtained by manual 
interactions with the display while the bottom two panels show the resulting amniotic fluid 
boundaries obtained from the instant invention's automatic application of 2D algorithms, 3D 
algorithms, combination heat and shock filter algorithms, and segmentation algorithms. 

5 After the contours on all the images have been delineated, the volume of the 

segmented structure is computed. Two specific techniques for doing so are disclosed in detail 
in U.S. Pat. No. 5,235,985 to McMorrow et al 9 herein incorporated by reference. This patent 
provides detailed explanations for non-invasively transmitting, receiving and processing 
ultrasound for calculating volumes of anatomical structures. 

10 Multiple Image Cone Acquisition and Image Processing Procedures: 

In some embodiments, multiple cones of data acquired at multiple anatomical sampling 
sites may be advantageous. For example, in some instances, the pregnant uterus may be too 
large to completely fit in one cone of data sampled from a single measurement or anatomical 
site of the patient (patient location). That is, the transceiver 10 is moved to different 
15 anatomical locations of the patient to obtain different 3D views of the uterus from each 
measurement or transceiver location. 

Obtaining multiple 3D views may be especially needed during the third trimester of 
pregnancy, or when twins or triplets are involved. In such cases, multiple data cones can be 
sampled from different anatomical sites at known intervals and then combined into a composite 
20 image mosaic to present a large uterus in one, continuous image. In order to make a composite 
image mosaic that is anatomically accurate without duplicating the anatomical regions 
mutually viewed by adjacent data cones, ordinarily it is advantageous to obtain images from 
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adjacent data cones and then register and subsequently fuse them together. In a preferred 
embodiment, to acquire and process multiple 3D data sets or images cones, at least two 3D 
image cones are generally preferred, with one image cone defined as fixed, and the other image 
cone defined as moving. 

5 The 3D image cones obtained from each anatomical site may be in the form of 3D 

arrays of 2D scanplanes, similar to the 3D array 240. Furthermore, the 3D image cone may be 
in the form of a wedge or a translational array of 2D scanplanes. Alternatively, the 3D image 
cone obtained from each anatomical site may be a 3D scancone of 3D-distributed scanlines, 
similar to the scancone 300. 

10 The term "registration" with reference to digital images means the determination of a 

geometrical transformation or mapping that aligns viewpoint pixels or voxels from one data 
cone sample of the object (in this embodiment, the uterus) with viewpoint pixels or voxels 
from another data cone sampled at a different location from the object. That is, registration 
involves mathematically determining and converting the coordinates of common regions of an 

15 object from one viewpoint to the coordinates of another viewpoint. After registration of at 
least two data cones to a common coordinate system, the registered data cone images are then 
fused together by combining the two registered data images by producing a reoriented version 
from the view of one of the registered data cones. That is, for example, a second data cone's 
view is merged into a first data cone's view by translating and rotating the pixels of the second 

20 data cone's pixels that are common with the pixels of the first data cone. Knowing how much 
to translate and rotate the second data cone's common pixels or voxels allows the pixels or 
voxels in common between both data cones to be superimposed into approximately the same x, 
y, z, spatial coordinates so as to accurately portray the object being imaged. The more precise 
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and accurate the pixel or voxel rotation and translation, the more precise and accurate is the 
common pixel or voxel superimposition or overlap between adjacent image cones. The precise 
and accurate overlap between the images assures the construction of an anatomically correct 
composite image mosaic substantially devoid of duplicated anatomical regions. 

5 To obtain the precise and accurate overlap of common pixels or voxels between the 

adjacent data cones, it is advantageous to utilize a geometrical transformation that substantially 
preserves most or all distances regarding line straightness, surface planarity, and angles 
between the lines as defined by the image pixels or voxels. That is, the preferred geometrical 
transformation that fosters obtaining an anatomically accurate mosaic image is a rigid 
10 transformation that doesn't permit the distortion or deforming of the geometrical parameters or 
coordinates between the pixels or voxels common to both image cones. 

The preferred rigid transformation first converts the polar coordinate scanplanes from 
adjacent image cones into in x, y, z Cartesian axes. After converting the scanplanes into the 
Cartesian system, a rigid transformation, T, is determined from the scanplanes of adjacent 
15 image cones having pixels in common. The transformation T is a combination of a three- 
dimensional translation vector expressed in Cartesian as t = (T x , T y , T 2 ), and a three- 
dimensional rotation R matrix expressed as a function of Euler angles fy, 0 Z around the x, y, 
and z axes. The transformation represents a shift and rotation conversion factor that aligns and 
overlaps common pixels from the scanplanes of the adjacent image cones. 

20 In the preferred embodiment of the present invention, the common pixels used for the 

purposes of establishing registration of three-dimensional images are the boundaries of the 
amniotic fluid regions as determined by the amniotic fluid segmentation algorithm described 
above. 
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Several different protocols may be used to collect and process multiple cones of data 
from more than one measurement site are described in FIGURES 13-14. 

FIGURE 13 illustrates a 4-quadrant supine procedure to acquire multiple image cones 
around the center point of uterine quadrants of a patient in a supine procedure. Here the 
5 patient lies supine (on her back) displacing most or all of the amniotic fluid towards the top. 
The uterus is divided into 4 quadrants defined by the umbilicus (the navel) and the linea-nigra 
(the vertical center line of the abdomen) and a single 3D scan is acquired at each quadrant. The 
4- quadrant supine protocol acquires four different 3D scans in a two dimensional grid, each 
corner of the grid being a quadrant midpoint. Four cones of data are acquired by the 
10 transceiver 10 along the midpoints of quadrant 1, quadrant 2, quadrant 3, and quadrant 4. Thus, 
one 3D data cone per uterine quadrant midpoint is acquired such that each quadrant midpoint is 
mutually substantially equally spaced from each other in a four-comer grid array. 

FIGURE 14 illustrates a multiple lateral line procedure to acquire multiple image cones 
in a linear array. Here the patent lies laterally (on her side), displacing most or all of the 
15 amniotic fluid towards the top. Four 3D images cones of data are acquired along a line of 
substantially equally space intervals. As illustrated, the transceiver 10 moves along the lateral 
line at position 1, position 2, position 3, and position 4. As illustrated in FIGURE 14, the inter- 
position distance or interval is approximately 6 cm. 

The preferred embodiment for making a composite image mosaic involves obtaining 
20 four multiple image cones where the transceiver 10 is placed at four measurement sites over the 
patient in a supine or lateral position such that at least a portion of the uterus is ultrasonically 
viewable at each measurement site. The first measurement site is originally defined as fixed, 
and the second site is defined as moving and placed at a first known inter-site distance relative 
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to the first site. The second site images are registered and fused to the first site images After 
fusing the second site images to the first site images, the third measurement site is defined as 
moving and placed at a second known inter-site distance relative to the fused second site now 
defined as fixed. The third site images are registered and fused to the second site images 
5 Similarly, after fusing the third site images to the second site images, the fourth measurement 
site is defined as moving and placed at a third known inter-site distance relative to the fused 
third site now defined as fixed. The fourth site images are registered and fused to the third site 
images 

The four measurement sites may be along a line or in an array. The array may include 
10 rectangles, squares, diamond patterns, or other shapes. Preferably, the patient is positioned 
such that the baby moves downward with gravity in the uterus and displaces the amniotic fluid 
upwards toward the measuring positions of the transceiver 10. 

The interval or distance between each measurement site is approximately equal, or may 
be unequal. For example in the lateral protocol, the second site is spaced approximately 6 cm 

15 from the first site, the third site is spaced approximately 6 cm from the second site, and the 
fourth site is spaced approximately 6 cm from the third site. The spacing for unequal intervals 
could be, for example, the second site is spaced approximately 4 cm from the first site, the third 
site is spaced approximately 8 cm from the second site, and the third is spaced approximately 6 
cm from the third site. The interval distance between measurement sites may be varied as long 

20 as there are mutually viewable regions of portions of the uterus between adjacent measurement 
sites. 

For uteruses not as large as requiring four measurement sites, two and three 
measurement sites may be sufficient for making a composite 3D image mosaic. For three 
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measurement sites, a triangular array is possible, with equal or unequal intervals. Furthermore, 
is the case when the second and third measurement sites have mutually viewable regions from 
the first measurement site, the second interval may be measured from the first measurement 
site instead of measuring from the second measurement site. 

5 For very large uteruses not fully captured by four measurement or anatomical sites, 

greater than four measurement sites may be used to make a composite 3D image mosaic 
provided that each measurement site is ultrasonically viewable for at least a portion of the 
uterus. For five measurement sites, a pentagon array is possible, with equal or unequal 
intervals. Similarly, for six measurement sites, a hexagon array is possible, with equal or 
10 unequal intervals between each measurement site. Other polygonal arrays are possible with 
increasing numbers of measurement sites. 

The geometrical relationship between each image cone must be ascertained so that 
overlapping regions can be identified between any two image cones to permit the combining of 
adjacent neighboring cones so that a single 3D mosaic composite image is produced from the 
1 5 4-quadrant or in-line laterally acquired images. 

The translational and rotational adjustments of each moving cone to conform with the 
voxels common to the stationary image cone is guided by an inputted initial transform that has 
the expected translational and rotational values. The distance separating the transceiver 10 
between image cone acquisitions predicts the expected translational and rotational values. For 
20 example, as shown in FIGURE 14, if 6 cm separates the image cones, then the expected 
translational and rotational values are proportionally estimated. For example, the (T x , T y , T z ) 
and 9y, 0 Z ) Cartesian and Euler angle terms fixed images p voxel values are defined 
respectively as (6cm, 0cm, 0cm) and (Odeg, Odeg, Odeg). 
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FIGURE 15 is a block diagram algorithm overview of the registration and correcting 
algorithms used in processing multiple image cone data sets. The algorithm overview 1000 
shows how the entire amniotic fluid volume measurement process occurs from the multiply 
acquired image cones. First, each of the input cones 1004 is segmented 1008 to detect all 
5 amniotic fluid regions. The segmentation 1008 step is substantially similar to steps 418-480 of 
FIGURE 7. Next, these segmented regions are used to align (register) the different cones into 
one common coordinate system using a Rigid Registration 1012 algorithm. Next, the registered 
datasets from each image cone are fused with each other using a Fuse Data 1016 algorithm to 
produce a composite 3D mosaic image. Thereafter, the total amniotic fluid volume is 
10 computed 1020 from the fused or composite 3D mosaic image. 

FIGURE 16 is a block diagram of the steps of the rigid registration algorithm 1012. 
The rigid algorithm 1012 is a 3D image registration algorithm and is a modification of the 
Iterated Closest Point (ICP) algorithm published by PJ Besl and ND McKay, in "A Method for 
Registration of 3-D Shapes," IEEE Trans. Pattern Analysis & Machine Intelligence, vol 14, 

15 no. 2, Feb. 1992, pp. 239—256. The steps of the rigid registration algorithm 1012 serves to 
correct for overlap between adjacent 3D scan cones acquired in either the 4-quadrant supine 
grid procedure or lateral line multi data cone acquisition procedures. The rigid algorithm 1012 
first processes the fixed image 1104 in polar coordinate terms to Cartesian coordinate terms 
using the 3D Scan Convert 1108 algorithm. Separately, the moving image 1124 is also 

20 converted to Cartesian coordinates using the 3D Scan Convert 1 128 algorithm. Next, the edges 
of the amniotic fluid regions on the fixed and moving images are determined and converted 
into point sets p and q respectively by a 3D edge detection process 1112 and 1132. Also, the 
fixed image point set, p, undergoes a 3D distance transform process 1116 which maps every 
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voxel in a 3D image to a number representing the distance to the closest edge point in p. Pre- 
computing this distance transform makes subsequent distance calculations and closest point 
determinations very efficient. 

Next, the known initial transform 1136, for example, (6, 0, 0) for the Cartesian T x , T y , 
5 T z terms and (0, 0, 0) for the^, fy, 0 Z Euler angle terms for an inter-transceiver interval of 6 
cm, is subsequently applied to the moving image by the Apply Transform 1140 step. This 
transformed image is then compared to the fixed image to examine for the quantitative 
occurrence of overlapping voxels. If the overlap is less than 20%, there are not enough 
common voxels available for registration and the initial transform is considered sufficient for 
10 fusing at step 1016. 

If the overlapping voxel sets by the initial transform exceed 20% of the fixed image p 
voxel sets, the q- voxels of the initial transform are subjected to an iterative sequence of rigid 
registration. 

A transformation T serves to register a first voxel point set p from the first image cone 
15 by merging or overlapping a second voxel point set q from a second image cone that is 
common to p of the first image cone. A point in the first voxel point set p may be defined as 
Pi =(*,-,,)>,•, z f .) and a point in the second voxel point set q may similarly be defined as 
<lj = ix J9 y J9 Zj) 9 If the first image cone is considered to be a fixed landmark, then the T 

factor is applied to align (translate and rotate) the moving voxel point set q onto the fixed voxel 
20 point set p. 

The precision of T is often affected by noise in the images that accordingly affects the 
precision of t and J?, and so the variability of each voxel point set will in turn affect the overall 
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variability of each matrix equation set for each point. The composite variability between the 
fixed voxel point set p and a corresponding moving voxel point set q is defined to have a cross- 
covariance matrix C M , more fully described in equation E8 as: 

C w =-E(A-P)(^-?) r E8 

5 where, n is the number of points in each point set and p and q are the central points in 

the two voxel point sets. How strong the correlation is between two sets data is determined by 
statistically analyzing the cross-covariance C pq . The preferred embodiment uses a statistical 
process known as the Single Value Decomposition (SVD) originally developed by Eckart and 
Young (G. Eckart and G. Young, 1936, The Approximation of One Matrix by Another of Lower 
10 Rank, Pychometrika 1,211-218). When numerical data is organized into matrix form, the SVD 
is applied to the matrix, and the resulting SVD values are determined to solve for the best 
fitting rotation transform R to be applied to the moving voxel point set q to align with the fixed 
voxel point set p to acquire optimum overlapping accuracy of the pixel or voxels common to 
the fixed and moving images. 

1 5 Equation E9 gives the SVD value of the cross-covariance C pg : 

C W = UDV' E9 

where D is a 3 x 3 diagonal matrix and U and V are orthogonal 3x3 matrices 

Equation E10 further defines the rotational R description of the transformation T in 
terms of U and V orthogonal 3x3 matrices as: 

20 R = UV T E10 
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Equation El 1 further defines the translation transform t description of the 
transformation T in terms of p ,q and R as: 

t = p-Rq Ell 

Equations E8 through Ell present a method to determine the rigid transformation 
5 between two point sets p and q - this process corresponds to step 1 152 in Figure 17. 

The steps of the registration algorithm are applied iteratively until convergence. The 
iterative sequence includes a Find Closest Points on Fixed Image 1 148 step, a Determine New 
Transform 1 152 step, a Calculate Distances 1 156 step, and Converged decision 1 160 step. 

In the Find Closest Points on Fixed Image 1 148 step, corresponding q points are found 
10 for each point in the fixed set p. Correspondence is defined by determining the closest edge 
point on q to the edge point of p. The distance transform image helps locate these closest 
points. Once p and closest -q pixels are identified, the Determine New Transform 1152 step 
calculates the rotation R via SVD analysis using equations E8-E10 and translation transform t 
via equation Ell. If, at decision step 1160, the change in the average closest point distance 
15 between two iterations is less than 5%, then the predicted-q pixel candidates are considered 
converged and suitable for receiving the transforms R and t to rigidly register the moving 
image Transform 1136 onto the common voxels p of the 3D Scan Converted 1108 image. At 
this point, the rigid registration process is complete as closest proximity between voxel or pixel 
sets has occurred between the fixed and moving images, and the process continues with fusion 
20 at step 1016. 

If, however, there is > 5% change between the predicted-q pixels and p pixels, another 
iteration cycle is applied via the Apply Transform 1140 to the Find Closest Points on Fixed 

1,1 llll I II I " 49 - Black Lowe&^graham 1 "* 

25315 DXUC-1-1020AP 816 Second Avenue 

x> a -nrKrr rJjnrv, Seattle, Washington 98104 

PATENT TRADEMARK OFFICE 206.381.3300 • F?206.381 .3301 



Image 1148 step, and is cycled through the converged 1160 decision block. Usually in 3 
cycles, though as many as 20 iterative cycles, are engaged until is the transformation T is 
considered converged. 

A representative example for the application of the preferred embodiment for the 
registration and fusion of a moving image onto a fixed image is shown in FIGURES 17A-17C. 

FIGURE 17A is a first measurement view of a fixed scanplane 1200A from a 3D data 
set measurement taken at a first site. A first pixel set p consistent for the dark pixels of AFV is 
shown in a region 1204A. The region 1204A has approximate x-y coordinates of (150, 120) 
that is closest to dark edge. 

FIGURE 17B is a second measurement view of a moving scanplane 1200B from a 3D 
data set measurement taken at a second site. A second pixel set q consistent for the dark pixels 
of AFV is shown in a region 1204B. The region 1204B has approximate x-y coordinates of 
(50, 125) that is closest to dark edge. 

FIGURE 17C is a composite image 1200C of the first (fixed) 1200A and second 
(moving) 1200B images in which common pixels 1204B at approximate coordinates (50,125) 
is aligned or overlapped with the voxels 1204A at approximate coordinates (150, 120). That is, 
the region 1204B pixel set q is linearly and rotational transformed consistent with the closest 
edge selection methodology as shown in FIGURES 13A and 13B from employing the 3D Edge 
Detection 1112 step.. The composite image 1200C is a mosaic image from scanplanes having 
approximately the same <j> and rotation Wangles. 

The registration and fusing of common pixel sets p and q from scanplanes having 
approximately the same <p and rotation 9 angles can be repeated for other scanplanes in each 
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3D data set taken at the first (fixed) and second(moving) anatomical sites. For example, if the 
composite image 1200C above was for scanplane #1, then the process may be repeated for the 
remaining scanplanes #2-24 or #2-48 or greater as needed to capture a completed uterine 
mosaic image. Thus an array similar to the 3D array 240 from FIGURE 5B is assembled, 
5 except this time the scanplane array is made of composite images, each composited image 
belonging to a scanplane having approximately the same <f> and rotation Wangles. 

If a third and a fourth 3D data sets are taken, the respective registration, fusing, and 
assembling into scanplane arrays of composited images is undertaken with the same 
procedures. In this case, the scanplane composite array similar to the 3D array 240 is 
1 0 composed of a greater mosaic number of registered and fused scanplane images. 

A representative example the fusing of two moving images onto a fixed image is shown 
in FIGURES 18A-18D. 

FIGURE 18A is a first view of a fixed scanplane 1220A. Region 1224A is identified as 
p voxels approximately at the coordinates (150, 70). 

15 FIGURE 18B is a second view of a first moving scanplane 1220B having some q 

voxels 1224B at x-y coordinates (300, 100) common with the first measurements p voxels at x- 
y coordinates (150, 70). Another set of voxels 1234A is shown roughly near the intersection of 
x-y coordinates (200, 125). As the transceiver 10 was moved only translationally, The 
scanplane 1220B from the second site has approximately the same tilt (f> and rotation Wangles 

20 of the fixed scanplane 1220A taken from the first lateral in-line site. 

FIGURE 18C is a third view of a moving scanplane 1220C. A region 1234B is 
identified as q voxels approximately at the x-y coordinates (250, 100) that are common with 
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the second views q voxels 1234 A. The scanplane 1220c from the third lateral in-line site has 
approximately the same tilt 0 and rotation 6 angles of the fixed scanplane 1220 A taken from 
the first lateral in-line site and the first moving scanplane 1220B taken from the second lateral 
in-line site. 

5 FIGURE 18D is a composite mosaic image 1220D of the first (fixed) 1220A image , 

the second (moving) 1220B image, and the third (moving) 1220C image representing the 
sequential alignment and fusing of q voxel sets 1224B to 1224A, and 1234B with 1234A. 

A fourth image similarly could be made to bring about a 4-image mosaic from 
scanplanes from a fourth 3D data set acquired from the transceiver 10 taking measurements at 
10 a fourth anatomical site where the fourth 3D data set is acquired with approximately the same 
tilt 0and rotation Wangles. 

The transceiver 10 is moved to different anatomical sites to collect 3D data sets by hand 
placement by an operator. Such hand placement could create the acquiring of 3D data sets 
under conditions in which the tilt $ and rotation 0 angles are not approximately equal, but 

1 5 differ enough to cause some measurement error requiring correction to use the rigid registration 
1012 algorithm. In the event where the 3D data sets between anatomical sites, either between a 
moving supine site in relation to its beginning fixed site, or between a moving lateral site with 
its beginning fixed site, cannot be acquired with the tilt (f> and rotation 6 angles being 
approximately the same, then the built-in accelerometer measures the changes in tilt (f> and 

20 rotation Wangles and compensates accordingly so that acquired moving images are presented if 
though they were acquired under approximately equal tilt ^ and rotation 0 angle conditions. 
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Demonstrations of the algorithmic manipulation of pixels of the present invention are 
provided in Appendix 1: Examples of Algorithmic Steps. Source code of the algorithms of 
the present invention is provided in Appendix 2: Matlab Source Code. 

While the preferred embodiment of the invention has been illustrated and described, as 
5 noted above, many changes can be made without departing from the spirit and scope of the 
invention. For example, other uses of the invention include determining the areas and volumes 
of the prostate, heart, bladder, and other organs and body regions of clinical interest. 
Accordingly, the scope of the invention is not limited by the disclosure of the preferred 
embodiment. 
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